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Abstract. DAMA is searching for rare processes by developing and using several kinds 
of radiopure scintillators: in particular, Nal(Tl), liquid Xenon and CaF2(Eu). Here only 
the results released so far on the WIMP annual modulation signature are summarized 
and compared with results from other experiments, including the recent re-analysis of 
CDMS-I data. Next perspectives are also shortly addressed. 

1 Introduction 

DAMA is devoted to the search for rare processes (such as WIMPs direct detec- 
tion, processes, charge-non-conserving processes, Pauli exclusion principle vi- 
olating processes, nucleon instabihty, solar axions and exotics |]l|,p|,pP,p|,p|,[7|P,p|, p^ , pdHT^ ) 
by developing and using low radioactive scintillators. The main experimental set- 
ups, which are running at present, are: the ~ 100 kg Nal(Tl) set-up, the ~ 6.5 
kg liquid Xenon (LXe) set-up and the so-called "R&D" apparatus. Moreover, a 
low-background germanium detector is operative underground since many years 
for measurements and selections of samples. 

In this paper only the results obtained in the search for WIMPs by exploiting 
the annual modulation signature with the ~ 100 kg Nal(Tl) set-up, their com- 
parison with those of other experiments (including the more recent re-analysis 
of the CDMS-I data) and some perspectives of the new LIBRA (Large sodium 
Iodine Bulk for RAre processes) set-up in preparation will be addressed. 

A full description of the ~ 100 kg Nal(Tl) set-up and of its performances can 
be found in ref. . It is worth to note that some upgrading has been performed 
since then; in particular, during August 2000 the electronic chain and the DAQ 
have been fully substituted achieving improved performances. 

The recoil/electron light ratio for ^^Na and ^^^I and the pulse shape discrim- 
ination capability have been measured by neutron source and upper limits on 
recoils have been measured in this set-up by exploiting the pulse shape discrim- 
ination technique H. Moreover, studies on possible diurnal variation of the low 
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energy rate in the data of the ~ 100 kg Nal(Tl) set-up have also been carried 
out. It could be expected because of the Earth's daily rotation; in fact, during 
the sidereal day the Earth shields a given detector with a variable thickness, 
eclipsing the "wind" of Dark Matter particles but only in case of high cross sec- 
tion candidates (to which small halo fraction would correspond). By analyzing a 
statistics of 14962 kg x day no evidence for diurnal rate variation with sidereal 
time has been observed this result supports that the effect pointed out by 
the studies on the WIMP annual modulation signature (see later) would account 
for a halo fraction > lO"'^ 

The main goal of the ~ 100 kg Nal(Tl) set-up is to investigate the WIMP 
annual modulation signature. In fact, the WIMPs are embedded in the galactic 
halo; thus, our solar system, which is moving with respect to the galactic system, 
is continuously hit by a WIMP "wind" which can be mainly searched for by 
WIMP elastic scattering on the target nuclei of the detector. In particular, since 
the Earth rotates around the Sun, which is moving with respect to the galactic 
system, it would be crossed by a larger WIMP flux in June (when its rotational 
velocity is summed to the one of the solar system with respect to the Galaxy) 
and by a smaller one in December (when the two velocities are subtracted). 
The fractional difference between the maximum and the minimum of the rate is 
expected to be of order of ~ 7%. 

The ~ 100 kg highly radiopure Nal(Tl) DAMA set-up || can effectively 
exploit such a signature because of its well known technology, of its high intrinsic 
radiopurity, of its mass, of its suitable control of all the operational parameters 
and of the deep underground experimental site. 

The annual modulation signature is very distinctive as we have already 
pointed out ||j5l|,@,|,|[|l^,|ll| . In fact, a WIMP-induced seasonal effect must 
simultaneously satisfy all the following requirements: the rate must contain a 
component modulated according to a cosine function (1) with one year period 
(2) and a phase that peaks around ~ 2"'' June (3); this modulation must be 
found in a well-defined low energy range, where WIMP induced recoils can be 
present (4); it must apply to those events in which just one detector of many 
actually "fires", since the WIMP multi-scattering probability is negligible (5); 
the modulation amplitude in the region of maximal sensitivity must be <7% (6). 
Only systematic effects able to fulfil these 6 requirements could fake this signa- 
ture; no one able to do that has been found or suggested on the contrary of 
what sometimes claimed by some author (see eg. rcf. Jl^). 

Results obtained by investigating the annual modulation signature in the 
data collected during the first four annual cycles have been released so far 
fHl^iR-iPl llol ^ lllll • The latter ones are summarized below and compared with those 
of some other experiments. At present the experiment is running taking data for 
the 7*'' cycle; at the end of this cycle the new LIBRA set-up (~ 250 kg of 
Nal(Tl)) wiU be installed. 
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A Model Independent Analysis of the 
Annual Modulation Data 



The annual modulation signature ofTers the possibility to obtain a model inde- 
pendent evidence for the presence of a WIMP component in the galactic halo. 
The large number of peculiarities (see above), which have to be satisfied for that, 
assures the absence of possible known systematic effects able to mimic such a 
signature (see above and ref. ^) and, therefore, an effective tool of investigation, 
as originally pointed out in ref. p[ . 

The data released so far refer to a total statistics of 57986 kg ■ day collected 
during four independent experiments of one year cycle each one. A model inde- 
pendent analysis of these data offers an immediate evidence of the presence of 
an annual modulation of the rate of the single hit events in the lowest energy 
interval (2-6 keV) as shown in Fig. |l|. There each data point has been obtained 
from the raw rate measured in the corresponding time interval, after subtract- 
ing the constant part. The results of the four years give consistent results with 
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Fig. 1. Model independent residual rate for single hit events, in the 2-6 keV cumulative 
energy interval, as a function of the time elapsed since January 1-st of the first year 
of data taking. The expected behaviour of a WIMP signal is a cosine function with 
minimum roughly at the dashed vertical lines and with maximum roughly at the dotted 
ones. 

proper period and phase. In particular, the test on the data of Fig. |^ disfavors 
the hypothesis of unmodulated behaviour (probability: 4- 10^^). No modulation 
is found in higher energy regions. 

We have extensively discussed the results of the investigations of all the pos- 
sible known systematics when releasing the data of each annual cycle; moreover, 
a dedicated paper |9j has been released on this topic. No known systematic ef- 
fect or side reaction able to mimic a WIMP induced effect has been found as 
discussed in details in ref. 0. In particular, no effect mentioned in ref. can 
be able to mimic it (see e.g. ref. [^). 

In conclusion, a WIMP contribution to the measured rate is candidate by 
the result of the model independent approach independently on the nature and 
coupling with ordinary matter of the possible WIMP particle. 
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3 Model Dependent Analyses of the 
Annual Modulation Data 

To investigate the nature and coupling with ordinary matter of a possible can- 
didate, a suitable energy and time correlation analysis is necessary as well as a 
complete model framework. We remark that a model framework is identified not 
only by general astrophysical, nuclear and particle physics assumptions, but also 
by the set of values used for all the parameters needed in the model itself and in 
related quantities (for example WIMP local velocity, vq, form factor parameters, 
etc.). 

For simplicity, initially we have considered the case of purely spin-independent 
coupled WIMP. In fact, often the spin- independent interaction with ordinary 
matter is assumed to be dominant since e.g. most of the used target-nuclei are 
practically not sensitive to SD interactions as on the contrary ^'^Na and ^^^I are 
and the theoretical calculations are even more complex. 

Moreover, the simplest model scenario has been considered as well as fixed 
parameters values Then, this case has been extended by considering the 

many uncertainties which exist on the astrophysical velocity distribution 
and the physical constraint which arises from the measured upper limits on re- 
coils Q . Then, some of the other possible particle scenarios have been considered 
such as extensions to the general case of WIMPs with both spin- independent (SI) 
and spin-dependent (SD) coupling and to that of WIMPs with inelastic scat- 
tering In these latter cases, the effect of the uncertainties on some other 
parameters has been included. Moreover, recently an investigation on the effect 
induced by different consistent halo models on the result for purely SI coupled 
WIMPs has also been carried out in ref. . 

At present the lightest supersymmetric particle named neutralino is con- 
sidered the best candidate for WIMP. Note, in particular, that the results of 
the data analyses |p|, p^[[TH summarized here and in the following hold for the 
neutralino, but are not restricted only to this candidate. 

In supersymmetric theories both the squark and the Higgs bosons exchanges 
give contribution to the coherent (SI) part of the neutralino cross section, while 
the squark and the Z'^ exchanges give contribution to the spin dependent (SD) 
one. Therefore, the differential energy distribution of the recoil nuclei in WIMP- 
nucleus elastic scattering can be calculated by means of the differential 

cross section of the WIMP-nucleus elastic processes: ^^^(^i^'-r) = (^^7) + 



y-j§j^j where v is the WIMP velocity in the laboratory frame and Ej^ is the 
recoil energy. 



3.1 WIMPs With Dominant SI Interaction in 
a Given Model Framev^rork 

As first scenario a full energy and time correlation analysis - properly account- 
ing for the physical constraint arising from the measured upper limit on recoils 
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[00 ~ been carried out in the framework of a given model for purely spin- 
independent coupled candidates with mass above 30 GeV. A standard maximum 
likelihood method has been used. Following the usual procedure we have built 
the y log-likelihood function, which depends on the experimental data and on 
the theoretical expectations for the considered model framework; then, y is min- 
imized and parameters' regions allowed at given confidence level are derived. 
Note that different model frameworks (see above) vary the expectations and, 
therefore, the cross section and mass values corresponding to the y minimum, 
that is also the allowed region at given C.L.. In particular, the inclusion of the 
uncertainties associated to the models and to every parameter in the models 
themselves as well as other possible scenarios largely enlarges the allowed region 
as discussed e.g. in ref. for the particular case of the astrophysical velocities. 
In the case considered in ref. the minimization procedure has been repeated 
by varying the WIMP local velocity, vq, from 170 km/s to 270 km/s to account 
for its present uncertainty. For example, the values mw = ij'^t.^) GeV and 
^(75/ = (5.7 ± 1.1) • 10"^ pb correspond to the position of y minimum when vq 
= 170 km/s, while mw = (431^^) GeV and £,asi = (5.4 ± 1.0) • lO'^ pb when 
vq = 220 km/s. Here, ^ is the WIMP local density in 0.3 GeV cm~^ unit, asi is 
the point-like SI WIMP-nucleon generalized cross section and mw is the WIMP 
mass. Fig. ^ shows the regions allowed at 3tT C.L. in such a model framework, 
when the uncertainty on is taken into account (solid contour) and when pos- 
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Fig. 2. A purely SI case: regions allowed at Scr C.L. on the plane S^crsi versus mw for a 
WIMP with dominant SI interaction and mass above 30 GeV in the model framework 
considered in ref. |^]: i) when vq uncertainty (170 km/s < vo < 270 km/s; continuous 
contour) has been included; ii) when also a possible bulk halo rotation as in ref. Q 
(dashed contour) is considered. Note that these regions hold for the given model frame- 
work (assumptions and parameters) and correspond to the superposition of all the 3 
a allowed regions obtained when varying vo in the allowed range; thus many "most 
likely" values correspond to it. Moreover, as widely known, the inclusion of present un- 
certainties on some other astrophysical, nuclear and particle physics parameters would 
enlarge these regions (varying again consequently the "most likely" values for each 
considered set). Some partial discussions can be found in [|7p|Jl^]. 
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sible bulk halo rotation is considered (dashed contour). For simplicity, no other 
uncertainty on the used parameters has been considered there (some of them 
have been included in the approach summarized in the next subsection pO||); 
thus, obviously the real allowed region including the effect of all the uncertain- 
ties is much larger than quoted for this simplified model framework as discussed 
at some extent also in ref. S and recently shown - as regards the halo model 
effect - in ref. |l| ( see Fig. ■ 

A quantitative comparison between the residuals and the modulation ampli- 
tudes for this particular model framework has been discussed in ref. More- 
over, we remark that the energy and time behaviour of the modulation ampli- 
tudes and residuals as well as the accounting for known uncertainties (such as 
quenching factors, form factors parameters, etc.) play a crucial role to obtain a 
meaningful comparison. 



10-8 _ 




10-11 ' ' ' ' ' — ' — ^ — ' — ^ — ' — ' — ' — ' — I 

bO 100 150 300 350 300 

(GeV) 

Fig. 3. Effect of halo modelling: the external continuous contour is the 3a annual 
modulation region for purely SI coupled WIMPs in the plane ^crsi (where asi = 
'^s'caiar""'' ) versus WIMP mass, obtained by considering all the (non-rotating) galactic 
halo models discussed in ref. [hd]. To point out the effect of the halo modelling alone, 
there is also shown the internal dashed region which corresponds to a simplified model 
of isothermal galactic halo in the same framework as ref. |8| , but assuming here the vq 
value at 220 km/s. 

In conclusion, the observed effect investigated in terms of a WIMP candidate 
with dominant SI interaction and mass above 30 GeV in the simplified model 
framework considered in ref. js), supports allowed WIMP masses up to about 
105 GeV (1 a C.L.) and even up to about 250 GeV (1 a C.L.) if different halo 
models are considered [|l2|. Lower £,(Tsi would be implied by the inclusion of 
known uncertainties on parameters (for example, on the quenching factors and 
on the form factor parameters) and on model features. 

Theoretical implications of these results in terms of a neutralino with domi- 
nant SI interaction and mass above 30 GeV have been discussed in ref. [T^Jl^ , 
while the case for an heavy neutrino of the fourth family has been considered 
in ref. ]l7[ . Let us comment that a correct evaluation and interpretation of the 
theoretical results are necessary; in fact, sometimes the uncertainties on the local 
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halo density are not taken into account and the rescahng procedures are not often 
apphed to correctly evaluate the WIMP local density. They must be considered 
for a reliable and correct presentation as well as the uncertainties which exist on 
some aspects of supersymmetric models in order to avoid that the given figure 
"would drive" the reader to wrong conclusion. Finally, note that the density of 
points in the calculated scatter plots do not represent a probability density. 



3.2 WIMPs With Mixed Coupling in Given Model Framework 

Since the ^'^Na and ^^^I nuclei are sensitive to both SI and SD couplings - on the 
contrary e.g. of ""'Ge and "'^*Si which are sensitive mainly to WIMPs with SI 
coupling (only 7.8 % is non-zero spin isotope in "°*Ge and only 4.7% of ^^Si in 
""'Si) - the analysis of the data has been extended considering the more general 
case 10 . This implies a WIMP having not only a spin- independent, but also 
a spin-dependent coupling different from zero, as it is also possible e.g. for the 
neutralino. 

Then, the log-likelihood function has been minimized - properly accounting 
also for the physical constraint set by the measured upper limit on recoils ||] - 
and parameters' regions allowed at given confidence level have been obtained. 
In particular, the calculation has been performed by minimizing the y function 
with respect to the £,asi, ^usd and mw parameters for each given 9 value. Here, 
asD is the point-like SD WIMP cross section on nucleon and tg9 is the ratio 
between the effective SD coupling constants on neutrons, Un, and on proton, a^; 
therefore, 9 can assume values between and tt depending on the SD coupling. In 
the present framework the uncertainties on vq have been included; moreover, the 
uncertainties on the nuclear radius and the nuclear surface thickness parameter 
in the SI form factor, on the b parameter in the used SD form factor and on the 
measured quenching factors |^ of these detectors have also been considered ||l^ . 
For simplicity. Fig. ^ shows slices for some mw of the region allowed at 3 u C.L. 
in the (^cts/i £,o'SD, "t-w) space for four particular couplings: i) 9 = (a„ =0 
and Op 7^ or \ap\ >> |a„|); ii) 9 = 7r/4 (op = a„); iii) 9 = Tr/2 (a„ ^ and Op 
= or |a„| >> \ap\); iv) 9 = 2.435 rad = -0.85, pure Z" coupling). The case 
Up = — a„ is nearly similar to the case iv). 

As already pointed out, when the SD contribution goes to zero (y axis in Fig. 

an interval not compatible with zero is obtained for £,asi- Similarly, when 
the SI contribution goes to zero (x axis in Fig. ||), finite values for the SD cross 
section are obtained. Large regions are allowed for mixed configurations also for 
^<^Si ^ 10~^ pb and £,asD < 1 pb; only in the particular case of 6* = -I (that is 
Op = and a„ 0) £,(Jsd can increase up to ~ 10 pb, since the ^'^Na and ^^^I 
nuclei have the proton as odd nucleon. Moreover, in ref. ||ic| we have also pointed 
out that: i) finite values can be allowed for £,asD even when £,asi — 3 • 10~^ 
pb as in the region allowed in the pure SI scenario considered in the previous 
subsection; ii) regions not compatible with zero in the ^ctsd versus mw plane are 
allowed even when ^asi values much lower than those allowed in the dominant 
SI scenario previously summarized are considered; iii) minima of the y function 
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Fig. 4. A mixed SI/SD case: example of slices of the region allowed at 3 (t C.L. in the 
{£,<^si, S.o'sD, rnw) space for some mw and 9 values in the model framework considered 
in ref. Only four particular couplings are reported here for simplicity: i) 9 = 0; ii) 
9 = 7r/4 iii) 9 = it/2; iv) 9 — 2.435 rad. Note that e.g. Ge experiments are sensitive 
mainly only to SI coupling and, therefore, cannot explore most of the DAMA allowed 
regions in this scenario. 

with both ^(7 SI a-nd ^asD different from zero are present for some mw and 9 
pairs; the related confidence level ranges between ~ 3 cr and ~ 4 cr pO[ |. 

Further investigations are in progress on these model dependent analyses 
to account for other known parameters uncertainties and for possible different 
model assumptions. As an example we recall that for the SD form factor an 
universal formulation is not possible since the internal degrees of the WIMP 
particle model (e.g. supersymmetry in case of neutralino) cannot be completely 
separated from the nuclear ones. In the calculations presented here we have 
adopted the SD form factors of ref. jlSj estimated by considering the Nijmengen 
nucleon-nucleon potential. Other formulations are possible for SD form factors 
and can be considered with evident implications on the obtained allowed regions. 

In conclusion, this analysis has shown that the DAMA data of the four annual 
cycles, analysed in terms of WIMP annual modulation signature, can also be 
compatible with a mixed scenario where both S,(7si and ^as d are different from 
zero. 



3.3 Inelastic Dark Matter 



It has been suggested |l9[| that the observed annual modulation effect could be 
induced by possible inelastic Dark Matter: relic particles that prefer to scatter 
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inelastically off of nuclei. The inelastic Dark Matter could arise from a mas- 
sive complex scalar split into two approximately degenerate real scalars or from 
a Dirac fermion split into two approximately degenerate Majorana fermions, 
namely x+ ^^^d x_, with a 5 mass splitting. In particular, a specific model fea- 
turing a real component of the sneutrino, in which the mass splitting naturally 
arises, has been given in ref. It has been shown that for the X- inelastic 
scattering on target nuclei a kinematical constraint exists which favours heavy 
nuclei (such as -'^^^I) with respect to lighter ones (such as e.g. "'°*Ge) as target- 
detectors media. In fact, X- can only inelastically scatter by transitioning to 
X+ (slightly heavier state than x-) ^^nd this process can occur only if the X- 

velocity is larger than vthr = \J mwN ^^^''C mwN is the WIMP-nucleus reduced 
mass (c — 1). This kinematical constraint becomes increasingly severe as the 
nucleus mass, mN, is decreased Moreover, this model scenario gives rise - 
with respect to the case of WIMP elastically scattering - to an enhanced mod- 
ulated component, 5'm, with respect to the unmodulated one, 5*0, and to largely 
different behaviours with energy for both Sq and Sm (both show a higher mean 
value) 

A dedicated energy and time correlation analysis of the DAMA annual modu- 
lation data has been carried out ||Tl[| handling aspects other than the interaction 
type as in ref. jl^ (in this way a particular model framework is fixed). 
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Fig. 5. An inelastic case: slices at fixed WIMP masses of the volume allowed at 3 o" 
C.L. in the space (^(Tp, 5, mw) obtained for the model framework considered in ref. 
PH]; some of the uncertainties on used parameters have been included |pd]]. Note that 
e.g. Ge experiments cannot explore most of the DAMA allowed regions in this scenario. 
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In this scenario of Dark Matter with inelastic scattering an allowed volume 
in the space {^ap,mw ,S) is obtained For simplicity, Fig. || shows slices of 
such an allowed volume at some given WIMP masses (3 a C.L.). It can be noted 
that when mw ^ nipf, the expected differential energy spectrum is trivially de- 
pendent on myy and in particular it is proportional to the ratio between ^cTp and 
mw', this particular case is summarized in the last plot of Fig. |^. The allowed 
regions have been obtained - as in the previous cases - by the superposition 
of those obtained when varying the values of the previously mentioned param- 
eters according to their uncertainties. Note that - as in the previous cases - 
each set of values (within those allowed by the associated uncertainties) for the 
previously mentioned parameters gives rise to a different expectation, thus to 
different "most likely" values. As an example we mention that when fixing the 
other parameters as in Ref. [|l^, the "most likely" values for a WIMP mass of 
70 GeV are: i) = 2.5 x 10"^ pb and (5 = 115 keV when vq = 170 km/s, it) 
^(jp = 6.3 X lO^** pb and 6 = 122 keV when vq = 220 km/s; they are in 5 re- 
gion were Ge and Si experiment are disfavoured. Finally, we note also here that 
significant enlargement of the given allowed regions should be expected when 
including complete effects of model (and related experimental and theoretical 
parameters) uncertainties. 

4 Proofs and Disproofs 

Let us preliminary remark that the claim for contradiction made by some au- 
thors (see e.g. pO| , ^ ) which use Ge target nuclei and different methodological 
approach is intrinsecally wrong. In fact (besides the usual uncertainties which 
always exist in the comparison of results achieved by different experiments) : i) 
the annual modulation signature gives a model independent evidence for WIMPs 
independently on their nature and coupling, while an exclusion plot is always 
model dependent; thus, no direct and model-independent comparison can be 
pursued among experiments which use different approaches, different techniques 
and, even more, different target nuclei; ii) within the same general assumptions 
for a model e.g. SI coupling and isothermal halo, the proper accounting for pa- 
rameters uncertainties, scaling laws uncertainties etc. significantly extends the 
allowed region and moves the best fit values; iii) there exist many model frame- 
works to which Na and I are sensitive and other nuclei (e.g. Ge and Si) are 
not. For example, a possible WIMP with a SI cross section of few 10^^ pb and 
with SD cross section of few 10"^ pb would produce a sizeable signal in DAMA 
but almost nothing in Ge and Si experiments. Moreover, also the given exclusion 
plots have a relative meaning since they are calculated under many assumptions; 
thus, they can drastically change either from model to model or just changing the 
value of each used experimental and theoretical parameter within its uncertainty. 

Let us shortly remind that some comments about other discussions of purely 
SD component as well as about the comparisons with other direct and indirect 
experiments can be found in ref. [p^ . Here, we only remind that the HEAT 
balloon experiment has confirmed an excess of high energy positrons in cosmic 
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rays which - if interpreted in terms of WIMP annihilation in a given particu- 
lar scenario - gives rise to a result compatible with that of DAMA within the 
uncertainties (see e.g. p2[). 

In the following only few comments on the CDMS-I result and on its com- 
parison with that of DAMA are given. 

As a general comment on the CDMS-I data as well as on their re-analysis 
in |2l], the arguments of our ref. |^ still hold. In addition, among the many cru- 
cial items necessary to credit the result itself and which are still missing, we cite: 

i) the measured energy spectra of gamma and electrons sources in running con- 
ditions. This can allow to show in a direct and safe way the detectors' response 
(from the peak positions and widths) and to derive the bulk detection efhciency; 

ii) the systematics associated to the used procedures since the measured counting 
rate is 60 counts keV'^ kg-^ d'^ (about 60 x 90keV x 15.8 kg x day = 85000 
events, for the considered selected data sample) and only few of them (23 in the 
re-analysis) survive the veto, the so-called "quality cuts" (made - as mentioned 
in ref. |2l| - to remove large periods with hardware troubles) and the discrimina- 
tion. In addition, these few events are mainly discarded by means of a neutron 
background modelling which suffers of some flaws, e.g.: a) the systematics due 
to the large amount of so-called "surface events" present in all the regions of 
the discrimination plot for the Si data (as shown by CDMS so far and now not 
included in Fig. 40 of [^); b) the uncertainties in the calculations based on the 
multiple scattering events. In fact, 4 events with both ionization yields Yl and 
Y2 < 0.5 (that is neutron candidates) and about 60-70 events with Y1,Y2 about 
1 (that is gamma/electron candidates) have been counted. The small number 
of gamma/electron multiple scattering has to be compared with the number of 
gamma/electron single scattering (after veto reduction procedure about 2 keV""'^ 
kg~-^ d~^ X 90keV x 15.8 kg x day — 2800 events). This is a relevant point 
since the neutron background modelling is based on the comparison of events 
with Y<0.5: 4 multiple-scatterings versus 23 single-scatterings. 

Furthermore, let us quickly comment about the comparison that generally is 
done - see e.g. Q - with the DAMA result. 

As mentioned above in many models Ge and Si target nuclei are insensitive 
to WIMP candidate to which instead Na and I are. Furthermore, in the models 
where all of them are sensitive, in order to obtain cross-section limits/regions 
from the experimental data, it is necessary to assume an astrophysical, particle 
and nuclear Physics model framework, to define a set of parameters values and 
to scale the cross sections on the different target nuclei to the one on proton. 
Therefore, every inferred result strongly depends e.g. on the chosen model, on 
the set of values fixed for the experimental and theoretical parameters and on 
the scaling laws for each involved nucleus. Thus, large uncertainty is associated 
to the results. This is a crucial point that must be accounted every time a com- 
parison among results from different experiments and with theoretical models is 
attempted, avoiding that a reader might have the wrong idea of the " universal- 
ity" of the conclusion. 
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In particular, as regards the purely spin-independent coupled WIMPs, which 
is practically the only scenario considered by other authors, the comparisons are 
often arbitrarily performed not with the region correctly endorsed by the DAMA 
collaboration (that is including the constraint from the measured upper limits) 
for the given simplified model framework of ref. ||] (see again |^). In fact, as 
an example, in Fig. ^ the exclusion limits claimed by CDMS and EDELWEISS 
are superimposed to the correct DAMA region for the simplified model with 
isothermal galactic halo, "simple" scaling laws and fixed assumptions for all the 
parameters as in ref. but assuming here the vq value at 220 km/s (that is 
neglecting the effect of its existing uncertainty on the allowed region, as correctly 
accounted in ref. ||] and here in the previous fig. ||). 




10° 10^ 10^ 10^ lO'* 
WTMPMass[GeV/c^] 



Fig. 6. A model dependent comparison for pur ely SI coupled WIMP: exclusion limits 
claimed by CDMS Q and EDELWEISS (24) , superimposed to the DAMA region 
allowed in the simplified model with isothermal galactic halo, "simple" scaling laws 
and fixed assumptions for all the parameters; in particular, here, for the comparison 
Vo has been set to the value 220 km/s (that the elTect of its existing uncertainty on 
the allowed region - correctly accounted in ref. ^ and here in the previous fig. ^ - has 
been neglected) . This figure is as in ref. [^l| , but here the correct region allowed by the 
cumulative DAMA data under the mentioned assumptions is shown. We further remind 
that the proper inclusion of the existing models' uncertainties significantly enlarges the 
allowed region (see e.g. the case of inclusion of vo uncertainties and of halo modelling 
in Fig. 1^ and ^ respectively). Also exclusion plots do not show "universal" boundary, 
but are affected by large experimental and theoretical uncertainties. 



Moreover, in ref. |21 a comparison is pursued by using some "most likely" 
values quoted in ref. however, it is worth to note that these "most likely" 
values have not a general validity, but depend on the specific assumptions (on 
halo, on particle, on nuclear physics) and experimental (e.g. quenching factors, 
detector response, etc.) and theoretical (form factors, scaling laws, etc.) param- 
eters' values. Varying them within the existing uncertainties varies the "most 
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likely" values. Realistic estimates of the expected number of counts in CDMS-I 
from the DAMA results for the purely SI coupled WIMPs can vary from about 
20 down to about 1. Moreover, let us remind that - even under the unrealistic as- 
sumption of no systematics in the hardware and software CDMS data reduction 

- the calculated exclusion plot has not an "universal" validity. 

A comparison between the Ge result and a particular model dependent ex- 
trapolation from the DAMA 2-6 keV cumulative residuals has also been ad- 
dressed in ref. |2^; this is again model dependent and, in particular, depends on 
the used quenching factors, form factor parameters, etc. In particular, the used 
procedures does take into account neither the energy behaviour of the measured 
modulation amplitude and the upper limits on recoils nor the proper values 
and uncertainties of parameters needed in the calculation. Thus, considerations 
similar to the previous case hold. 

Finally, in Fig. 47 of ref. |^ - as sometimes happens in other presentations 

- the theoretical expectations for the particular case of the neutralino are shown 
without taking into account the uncertainties on the local halo density and the 
rescaling procedures (as well as on some other astrophysical, nuclear and particle 
physics aspects). This incorrect procedure leads to wrong conclusion. 

Summarizing, there is no solid scientific reason for the CDMS claim for con- 
tradiction whatever scenario would be considered. 



5 Towards LIBRA 

At present our main efforts are devoted to the realization of LIBRA (Large 
sodium Iodine Bulk for RAre processes in the DAMA experiment) consisting of 
~ 250 kg of Nal(Tl). New radiopure detectors by chemical/physical purification 
of Nal and Tl powders as a result of a dedicated R&D with Crismatec have been 
realized. The installation of this set-up is foreseen in fall 2002; this will allow 
us to increase the sensitivity of the experiment. Some related arguments can be 
found in ref. [p5| . 

6 Conclusions 

The DAMA annual modulation data of four annual cycles |4||^,|7|,p|j9| JTo| , [ri| have 
been analysed by energy and time correlation analysis in terms of purely SI, 
purely SD, mixed SI/SD, "preferred" inelastic WIMP scattering model frame- 
works. 

To effectively discriminate among the different possible scenarios further in- 
vestigations are in progress. In particular, the data of the 5*'* and 6*^ annual 
cycles are at hand, while the set-up is running to collect the data of a 7*'* an- 
nual cycle. Moreover, the LIBRA (Large sodium Iodine Bulk for RAre processes) 
set-up is under construction to increase the experimental sensitivity. 
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